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ABSTRACT 


An  infinitely  long,  ring-stiffened,  submerged,  elastic  cylinder 
having  uniformly  spaced  elastic  bulkheads  is  the  structure  considered. 
Loading  is  applied  by  a plane  acoustic  shock  wave  with  front  parallel 
to  the  cylinder  axis.  Dynamic  pressure  in  the  fluid  is  resolved  into 
a free-field  incident  part  and  a scattered  part.  Structural  response 
emd  scattered  pressure  in  the  surrounding  fluid  are  found  using  finite 
element  modeling  of  structure  and  fluid.  Introduction  of  Fourier 
series  makes  the  fluid  region  mathematically  two-dimensional.  A radia- 
tion, or  non-reflecting,  condition  at  the  outer  boundary  of  the  fluid 
region  is  shown  to  give  good  results.  A p2u:ametric  study  is  made  of 
effects  of  shock  pulse  rise  time  and  duration  on  structural  response. 
Results  are  presented  as  combinations  of  shock  pressure  and  submergence 
pressure  just  sufficient  to  induce  structural  failure. 


! 

f 


TABLE  OF  CONTENTS 


I.  INTBODUCTKW 8 

II.  WAVE  PROPAGATION 10 

A.  PRESSURE  RESOLUTION 10 

B.  RADIATION  BOUNDARY  CONDITION 13 

C.  ONE-DIMEICIONAL  ANALYSIS 14 

D.  TWO-DIMENSICMJAL  ANALYSIS 15 

1.  Boundeury  Conditions  - --  --  --  --  --  --  --  -15 

2.  Fourier  Analysis  - --  --  --  --  --  --  --  --  -16 

III.  STRUCTURAL  RESPWSE 20 

A.  STRUCTURAL  MODEL 20 

B.  FLUID  MODEL 21 

C.  SOLUTION  TECHNIQUE 22 

D.  MODEL  PARAMETERS  23 

E.  FAILURE  CRITERIA 23 

IV.  RESULTS 24 

A.  INTRODUCTION  . 24 

B.  EFFECTIVENESS  OF  RADIATION  BOUNDARY  24 

C.  INFLUENCE  OF  PULSE  RISE  TIME  AND 

DURATION  ON  STRUCTURAL  FAILURE  26 

V.  CONCLUSIONS  40 

APPENDIX  A.  FLUID  COEFFICIENT  MATRICES  <1 

APPENDIX  B.  NUMERICAL  RESULTS  FOR  FAILURE  BOUNDARIES  44 

LIST  OF  REFERENCES -___-57 

INITIAL  DISTRIBUTION  LIST 58 


J 


5 


LIST  OF  FIGURES 


1.  Wave  front  - 

2.  Pressiore-time  history  of  a cubic  ramp-box  wave 

3.  Pareuneters  used  in  Fourier  analysis  ------ 

4.  Region  represented  by  finite  element  model  - - 
Effective  pulse  duration  ----------- 

Failure  boundaries  for 

Failure  boundaries  for 

Failure  boundfuries  for 

Failure  boundaries  for 

Failure  boundaries  for 

Failure  boundaries  for 

Failure  boundaries  for 

Failure  boundaries  for  — ^ * 0.8,  = 4 

a a 

cT 

Failure  boundaries  for 


5. 

6. 

7. 

8. 

9. 

10. 

11. 

12. 

13. 

14. 

15. 

16. 


n 1 

cT 

a 

U • X f 

a 

CT 

r 

cT 

a 

U*  A f 

a 

CT 

r 

a 

0.4, 

cT 

a 

n Q 

cT 

a 

U • O f 

a 

cT 

r 

n 1 

cT 

a 

w • X f 

a 

cT 

r 

cT 

a 

U • A ^ 

a 

cT 

r 

n A 

cT 

a 

a 

cT 

r 

0.8, 

cT 

a 

a 

cT 

r 

cT 

a 

U • X ^ 

a 

12 

17 

17 

20 

26 

29 

30 

31 

32 

33 

34 

35 

36 

37 


CT 


Failure  boundauries  for 


cT 


a 

cT 


0.2, 0.4, 0.8,  — -8 38 

A 


Failure  boundaries  for  — — - 0.97,  — ■ » 

A A 


39 


ACKNOWLEDGEMENT 


The  author  wishes  to  express  his  deep  appreciation  to  Dr.  Robert 
E.  Newton,  Professor  of  Mechanical  Engineering,  for  his  guidemce  and 
help  throughout  this  investigation.  A special  note  of  gratitude  is 
extended  to  the  staff  of  the  W.  R.  Church  Computer  Center  of  the  Naval 
Postgraduate  School. 

Finally,  the  author  thanks  his  feuoily.  Ana  Maria,  Jose  Gustavo, 
Maria  Alejandra  and  Ricardo  Alfonso  for  their  love  and  patience. 


.1 

I 


I.  INTRODUCTION 


nils  investigation  is  the  third  in  a sequence  of  finite  element 
studies  of  fluid-structure  interaction.  The  initial  study  was  by 
Atchison  [1]  and  the  second  by  Newton  and  Atchison  [2].  Tlie  system  c(xi- 
sidered,  as  in  the  earlier  studies,  is  a submarine  subjected  to  a shock 
wave  resulting  from  an  underwater  explosion.  The  submarine  is  modeled 
as  a ring-stiffened  cylinder  of  infinite  length  with  equally  spaced 
elastic  bulkheads.  The  shock  wave  is  a plane  acoustic  wave  with  front 
parallel  to  the  longitudinal  axis  of  the  cylinder.  The  structure- fluid 
equations  are  coupled  by  terms  resulting  from  dynamic  interaction  at 
the  interface  boundary. 


A.  PRIOR  STUDIES 

It  is  shown  by  both  authors  in  references  [1]  and  [2]  that  the  three- 
dimensional  problem  may  be  resolved  into  two-dimensional  subproblems, 
acoustic  wave  propagation  and  transient  response  of  a submerged  struc- 
ture, using  the  finite  element  method. 

Acoustic  wave  propagation  and  spatial  fluid  discretizaticm  are  pre- 
sented in  [1] . The  superposition  theorem,  which  permits  the  rigid  body 
pressure  distribution  to  be  obtained  from  a two-dimensional  analysis  is 
also  presented  in  [1] . 

A refined  structural  model  is  used  in  [2].  m this  paper,  the  cir- 
cumferential variations  of  fluid  pressures  auid  structure  displacements 
are  represented  by  Fourier  expansions.  Fluid  pressures,  structural  de- 
flections, shell  and  ring  stresses  are  tabulated  for  each  harmonic 
through  n>4.  Extreme  stresses  resulting  from  superposition  of  harmonic 
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contributions  are  given  also.  Two  separate  failure  criteria  are  used 


to  determine  critical  ccmtbinations  of  submergence  pressure  and  shock  ^ 
pressure.  The  first  of  these  utilizes  limit  analysis  to  evaluate  combined  ] 
effects  of  longitudinal  and  eucial  stress  in  the  shell  wall  at  a bulk-  I 
head.  The  second  postulates  that  shell  collapse  will  occur  when  the  | 
von  Mises  yield  condition  is  satisfied  away  from  the  bulkhead. 


B.  PRESENT  STUDY 

In  the  present  study  the  need  for  a separate  solution  of  the  acous- 
tic propagation  is  eliminated,  thus  effecting  a significant  reduction 
in  computation.  This  improvement  is  made  by  dealing  sepaurately  with 
the  incident  shock  pressure  - am  entity  which  is  unaffected  by  the  pres- 
ence of  the  structure.  The  remainder  of  the  dynamic  fluid  pressure 
(called  the  "scattered"  pressure)  is  represented  in  the  finite  element 
modeling  of  fluid-structure  response.  As  in  [2],  the  azimuthal  vauria- 
tion  of  fluid  pressure  and  structural  displacement  are  represented  by 
trigonometric  series  so  that  the  fluid  region  modeled  is  mathematically 
two-dimens ional . 

A further  improvement  introduced  here  is  the  imposition  of  a radia- 
tion boundary  condition  at  the  cylindric  outer  boundary  of  the  fluid. 
Although  this  condition  is  necessarily  approximate,  it  is  found  to  re- 
duce significantly  the  effects  of  boundary  reflections  on  structural 
response. 

A systematic  study  is  mauie  of  the  effects  of  shock  pulse  duration 
and  rise  time  on  structural  response.  Results  are  presented  in  the  form 
of  critical  ccmbinations  of  shock  pressure  and  submergence  pressure  just 
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I I . WAVE  PROPAGATION 


A.  PRE*?URE  RESOLUTION 

A s.,  lerged  structure  subjected  to  an  underwater  shock  wave  experi- 
ences, at  a given  time,  a total  dynamic  pressure  which  can  be  considered 
to  have  three  components 


P = Pi  + Pj^  + Pj.  (1) 

where 

p is  the  total  dynamic  pressure  in  the  fluid , 
p^  is  the  pressure  due  to  the  incident  wave,  considering 
the  structure  absent , 

Pj^  is  the  pressure  due  to  the  reflected  wave,  consider- 
ing no  structural  deflection, 

Pj,  is  the  pressure  radiated  by  the  structure  due  to 
structural  deflection. 

There  are  two  useful  combinations  of  these  pressures  for  purposes 
of  analysis.  Consider  first 

Prig  ” Pi  ^ Pr  (2) 

where  Pj-ig  represents  the  total  dynamic  pressure  field  with  a rigid 
structure.  Next  consider 

Ps  “ Pr  ^ Pr 

where  p^  represents  the  scattered  pressure,  which  includes  both  re- 
flected and  radiated  parts. 


Instead  of  considering  the  total  dynamic  pressure  to  be  the  sum  of 

p . and  p as  in  [1]  and  [2] , the  resolution 
rig  r 


p = p.  + p 


will  be  used.  This  resolution  has  been  employed  widely  in  similar  prob- 
1ms,  e.g.,  by  Mindlin  and  Bleich  [3],  Geers  [4]  and  Everstine  [5]. 

The  discretized  equations  for  the  structure  and  fluid  given  in  [1] 
and  [2]  are 


M^+K^*L£ 


££+D£+H£=-PL  6 


where  6 is  a column  vector  of  structural  nodal  displacements  and  £ 
is  a column  vector  of  fluid  nodal  pressures.  M,  K,  D and  li  are  sym- 
metric coefficient  matrices,  L is  a matrix  describing  coupling  at  the 
structure-fluid  interface,  and  p is  the  fluid  density.  Superior  dots 
denote  time  derivatives  and  superscript  T denotes  tremsposition. 

Considering  first  the  incident  pressure  only,  a time-dependent  column 
vector  ^ is  added  to  correct  for  the  unsatisfied  boundary  condition 
at  the  fluid-structure  interface.  It  is  required  that 


2 Ei  + £ ii  + H Ei  =•  1 


Considering  next  the  scattered  pressure,  it  is  required  that 
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The  equation  for  the  structure  (5)  becomes 


P 


I 


m6  + k6  = L£^  + L£^  (9) 

Therefore,  equations  (7) , (8)  £Uid  (9)  are  the  new  equations  for  the  con- 
plete  elasto-hydrodynamic  problmn  formulated  by  Atchison  [1]  and  Newton 
and  Atchison  [2] . 

Rather  than  use  equation  (7)  to  find  the  relation  between  the  in- 
cident pressure  and  the  fluid  particle  velocity  normal  to  the  boundary 
may  be  used. 

Consider  a plcine  wave  traveling  in  the  negative  x-direction.  The 
free-field  incident  pressure  is  given  by 

p^  = g(x  + ct)  (10) 

where  c is  the  acoustic  velocity  of  propagation.  Fig.  1 represents 
the  wave  front. 


Pi 


Applying  Newton's  second  law  to  a fluid  particle,  it  can  be  shown 
that  the  fluid  acceleration  is  given  by 


1 ^^i  _ _ £' 


p dx 


(11) 
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However 


Pi  =-5t- 


= g'c 


Therefore,  in  (11) 


1 • 

P. 

pc 


(12) 


Formula  (12)  gives  the  relation  between  the  incident  pressure  and  the 
fluid  particle  velocity.  The  complete  problem  now  is  reduced  to  solv- 
ing the  set  of  linear  equations  given  by  (8)  and  (9)  using  the  result 
est^U^lished  by  equation  (12) . 

Fourier  representations  are  employed  in  deriving  the  discretized 
counteiqjarts  of  equations  (8)  and  (9)  for  each  pressure  harmonic  and 
displacement  heunnonic.  Because  both  the  structure  and  the  fluid  region 
are  axisymmetric,  there  is  no  interharmonic  coupling. 

Therefore,  for  the  structure  displacement,  let 


6 =6  cos  n0 

~ ~n 


(13) 


and,  for  the  fluid  pressure,  let 


£ = cos  n0  (14) 

where  n-0,1,2,...  and  summation  with  respect  to  n is  implied. 

For  each  Fourier  order  n there  is  a separate  pair  of  equations  (8) 

and  (9)  governing  the  corresponding  pair  5 and  p of  displacement 

~n  “"n 

and  pressure  vectors.  The  coefficient  matrices  of  these  equations  are 
order-dependent.  Formulas  for  these  matrices  are  given  in  Appendix  A. 


B.  RADIATION  BOUNDARY  CONDITION 

Finite  element  modeling  of  infinite  regions  always  requires  seme 

strategem  to  avoid  significant  errors  due  to  substitution  of  a finite 
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region.  In  most  static  problems  and  some  steady-state  dynamic  problems 
it  suffices  to  place  the  boundaries  "far  enough"  frcan  the  region  of 
interest.  This  scheme  is  also  appliccdsle  to  the  present  problem  featur- 
ing wave  propagation  in  the  fluid.  Here  the  requirement  of  "far  enough" 
is  satisfied  if  waves  traveling  from  structure  to  outer  boundeury  and 
back  to  the  structure  airrive  after  the  occurrence  of  m^ucimum  stresses  in 
the  structure.  Meeting  this  requirement  led,  in  the  present  applica- 
tion, to  a finite  element  model  in  which  over  80  percent  of  the  degrees 
of  freedom  represent  fluid  pressures. 

An  alternate  means  for  simulating  an  infinite  region  utilizes  a 
radiation,  or  non- reflecting,  condition  at  the  outer  boundeuries.  Such 
a condition  was  first  proposed  by  Zienkiewicz  and  Newton  [6]  in  a form 
which  is  exact  for  plane  waves  normally  incident  upon  the  boiuidary. 
Chenault  successfully  ^ployed  the  condition  in  a study  of  two-dimensional 
added  mass  and  damping  [7] . Dean  constructed  and  successfully  demonstra- 
ted an  extension  to  spherically  spreading  waves  [8] . For  both  plane 
euid  spherical  waves  the  corresponding  non-reflecting  conditions  are 
exact,  but  no  such  result  is  possible  for  two-dimensional  (cylindric) 
spreading.  Despite  the  necessity  for  approximation,  Bai  constructed  and 
used  a radiation  condition  for  two  dimensions  [9].  An  extension  of  Bai's 
scheme  is  developed  in  the  section  entitled  "Two-Dimensional  Analysis". 

C.  ONE-DIMENSIONAL  ANALYSIS 

A one-dimensional  analysis  using  the  new  pressure  resolution  has  been 
done  in  order  to  compare  new  results  with  those  previously  obtained  by 
Atchison  [1] . 

The  region  considered  is  a semi-infinite  fluid  strip  of  unit  cross 
sectional  area.  Through  this  region  a shock  wave  of  known  shape  is 
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propagating.  The  semi-infinite  strip  is  modeled  with  a one  degree  of 
freedom  structure  at  one  end  and  the  radiation  boundary  condition  is 
applied  at  the  other  end. 

The  results  obtained  were  compcured  with  those  obtained  by  using 
CSMP  (continuous  system  modeling  program) . Good  agreement  was  found. 
Accordingly,  the  Houbolt^ integration  method  [10],  the  line2u:  element, and 
the  restrictions  for  node  spacing  and  time  step  as  used  in  [1]  and  [2] 
were  found  to  give  excellent  results  when  used  with  the  new  pressure 


resolution. 


0.  TWO-DIMENSIONAL  ANALYSIS 


1.  Boundary  Conditions 


A non-reflecting  or  radiation  boundary  condition  is  presented 
and  used  with  this  fluid- structure  investigation. 

If  it  is  wished  to  let  a plane  wave  traveling  in  the  positive  x- 
direction  pass  out  of  the  region  being  investigated,  the  condition  re- 
quired is 


3p  1 • 

If  - - - P (15) 

In  dealing  with  two-dimensional  problems,  the  axisymmetric  ap- 
proximation proposed  by  Bai  [9]  is  adapted.  The  (scattered)  pressure 
distribution  for  a wave  moving  in  the  positive  r-direction  is  assumed  to 


be  described  by 


■ (r- 


g (r-ct)cos  n8 
n 


where,  as  before,  there  is  summation  with  respect  to  index  n. 
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The  boundary  condition  imposes  a requirement  on  the  normal  deri- 
vative of  the  pressure  at  the  boundary.  The  partial  derivative  is  cal- 
culated from  equation  (16)  as 

1 —3/2  —1/2 

“ “ 2 g^(r-ct)  cos  nO  + r"  ' g^r-ct)  cos  n0  (17) 


Also 


. 3Ps 

•s  -sr 


■s  -1/2 

P„  * ^ 9'  (i^”Ct)  cos  n6(-c) 


or 


c 


-1/2  , , 

- r g'(r- 


g'(r-ct)  cos  n0 
n 


Therefore,  equation  (17)  may  be  rewritten  as 


^Ps 

jr 


s 

2r 


s 

c 


(18) 


This  formula  will  be  introduced  as  the  radiation  boundary  condition  for 
the  solution  of  the  interaction  problem.  Data  concerning  the  usefulness 
of  this  approximation  are  given  in  the  section  titled  Results. 

2.  Fourier  Analysis 

It  is  required,  in  order  to  solve  equations  (8)  and  (9),  to  know 
both  the  incident  pressure  auid  the  normal  component  of  the  fluid  acce- 
leration at  the  structure-fluid  interface. 

The  pressure-time  history  of  the  cubic  ranp-box  wave  used  may  be 

represented  as  shown  in  Fig.  2,  where  T is  the  pulse  rise  time,  x 

r d 

is  the  pulse  duration  and  i is  the  subsidence  time. 

For  a given  incident  pressure  wave,  characterized  by  rise  time, 
pulse  duration,  subsidence  time,  acoustic  velocity  and  fluid  density,  it 
is  required  to  find  Fourier  con^jonents  for  the  incident  pressure  and 
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Fig.  2.  pressure-time  history  of  a cubic  ramp-box  wave. 


fluid  acceleration  at  the  shell  surface  using  a chosen  time  step  h. 
These  Fourier  components  are  needed  for  as  many  time  steps  as  are  used 
in  solving  the  fluid-structure  equations. 

Basic  parameters  for  the  Fourier  analysis  cure  shown  in  Fig.  3. 


Fig.  3.  Parameters  used  in  Fourier  analysis 


For  the  cubic  wave  front,  pressure  and  pressure  derivative  fonnul^ls  are 
derived  by  Atchison  in  [1] 


2 2 * 
OC  - 2C  )p  , 


Pi 


(C  - C ) , 


t < T 


T < t < (T  + T.) 
r — — r d 


t < T 


T < t < (T  + T,) 
r — — r d 


(19) 


(20) 


where  C * — and  p*  is  the  shock  pressure, 
r 

The  subsidence  is  described  by  the  same  kind  of  cubic  transition 
as  the  rise. 

Relations  for  the  normal  (perpendicular)  component  v^  of  the 
fluid  acceleration  can  be  derived  using  equations  (12)  and  (19) . Thus 


(C  - COS0  , t < 


(21) 


, Tr  < t < (Tr  + T^) 


Therefore,  equations  (19)  and  (21)  are  used  as  basic  formulation 
for  the  Fourier  analysis. 

Fourier  conqponents  for  the  incident  pressure  and  fluid  accelera- 
tion are  calculated  for  each  time  step  using  24  intervals  over  the 
cylinder  from  6 > O'*  to  6 > 180*.  It  is  not  necesseury  to  include  values 
of  9 from  180*  to  360*  because  of  symmetry  about  6 > 0*. 
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where  w ■ w„  - 1/2  emd  all  other  w « 1. 

ON  J 

For  the  fluid  normal  acceleration  con^onents 


N 

J=0  ^ 

N 

®n  “ n9j 

J=0 


N - 24  (23) 

n - 1,2, 3, 4 


Required  values  of 


are  found  from  equations  (19) 


and  (21) . 


New  Fourier  resolutions  are  required  when  any  characteristic 


time  of  the  shock  pulse  is  changed. 
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III.  STRUCTURAL  RESPONSE 

A.  STRUCTURAL  MODEL 

A refined  model  is  used  for  the  structure.  A total  of  eight  struc- 
tural elements  and  27  structural  degrees  of  freedom  per  Fourier  component 
represent  the  final  mesh  and  it  is  identical  with  the  one  utilized  by 
Newton  and  Atchison  [2] . 

As  described  in  [2] , the  structure  is  a ring-stiffened  elastic  shell 
of  infinite  length  with  elastic  bulkheads  at  spacing  2L.  Fig.  4 shows 
a cross-section  of  the  repeating  length  of  structure  and  the  associated 
fluid  region. 


zi 


r 


Fig.  4.  Region  represented  by  finite  element  model. 

It  is  assumed  that  the  plcuie  of  a bulkhead,  z»0,  and  the  plane  mid- 
way between  bulkheads,  z>«L,  eu:e  planes  of  synmetry  for  the  structure 
and  the  response  displacmnents . 

As  in  reference  [2] , the  effects  of  the  stiffening  rings  are  in- 
cluded by  treating  the  shell  as  orthotropic.  Stress  calculations  take 
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into  account  the  geometry  of  the  rings.  The  effective  shell  thickness 


in  the  circumferential  direction  includes  a contribution  by  the  rings. 

The  bulkheads  furnish  elastic  restraint  against  radial  displacement 
of  the  shell.  This  is  accomplished  by  treating  the  bulkhead  as  a 
diaphragm  of  uniform  thickness. 

As  developed  in  reference  [2]  a geometric  stiffness  matrix  modifies 
the  structural  stiffness.  It  considers  the  effects  of  the  membrane 

compressive  stresses  of  the  axisymmetric  mode  acting  to  amplify  deflec- 
tions in  the  higher  harmonics  (n>l) . 

Further  information  concerning  the  structural  model  is  given  in 
reference  [2] . Details  concerning  application  of  Fourier  series  to 
analysis  of  eucisymmetric  shells  are  given  by  Grafton  atnd  Strome  [11] . 

B.  FLUID  MODEL 

As  indicated  in  Fig.  4 the  length  of  the  fluid  region  is  one-half 

of  the  bulkhead  spacing,  its  inner  radius  is  the  shell  radius  "a",  and 

the  outer  radius  is  r . 

o 

Rectangular  elements  having  corner  nodes  and  linear  shape  function 
are  used.  Along  the  fluid-structure  interface  a fluid  node  is  located 
at  each  structural  node. 

The  radiation  boundary  condition  as  developed  previously  is  applied 
at  the  outer  radius  r^.  For  the  purpose  of  studying  the  effectiveness 
of  the  radiation  boundzury  condition,  the  same  combination  of  pulse  para- 
meters and  submergence  pressure  has  been  employed  using  a>5  meters  and 

a sequence  of  values  of  r from  8.15  meters  to  17.6  meters. 

o 

In  finding  critical  combinations  of  shock  amd  subnergence  pressure 
to  induce  structural  failure  two  different  meshes  are  utilized.  In  each 


4- 


fl 


1 i 


1 
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case  the  fluid  region  has  inner  radius  a = 5 meters  and  outer  radius  r 

o 

16.025  meters.  The  first  mesh  has  90  fluid  elements  and  120  fluid 
nodes  and  each  elment  has  a thickness  (radial)  of  0.525  meters.  The 
second  has  174  fluid  elements  and  225  fluid  nodes  with  element  thickness 
0.2625  meters.  The  latter  mesh  is  used  when  the  rise  time  of  the  wave 
front  is  small  in  order  to  meet  the  restriction  given  by  Atchison  [1]  as 


where  s is  the  node  spacing.  Two  different  time  steps  are  employed 
depending  on  the  size  of  the  mesh.  For  the  finer  mesh  a tire  step  of 
0.1  milliseconds  is  used.  For  the  coarser  mesh  the  time  step  is  0.3 
milliseconds. 


C.  SOLUTION  TECHNIQUE 

The  solution  process  is  very  similar  to  the  one  used  by  Newton  and 
Atchison  [2]. 

There  is  a sepaurate  solution  for  each  of  the  five  harmonics  considered 
(from  n*0  to  n=4) . Once  the  solution  for  the  first  harmonic  is  completed, 
the  solution  for  the  second  starts  and  the  other  harmonics  follow  in 


order. 


Fourier  components  for  the  incident  pressure  and  normal  fluid 


acceleration  are  introduced  as  input  for  each  solution. 

For  any  harmonic  the  procedure  employed  first  advances  by  one  time 
step  using  equation  (9) . The  value  of  at  the  forward  step  is  obtained 
by  parabolic  extrapolation.  With  ^ and  ^ at  the  forward  point,  equa- 
tion (8)  is  then  used  to  meUce  the  forward  step  for  - Successive 
iterations  are  performed  using  the  latest  values  of  £g  and  ^ on  the  right 
hand  sides  of  equations  (9)  emd  (8) . 
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Detailed  information  eUaout  this  solution  technique  is  given  in  refer- 


ence [2] . 


D.  MODEL  PARAMETERS 


The  following  parameters  are  used  for  structural  response  calcula- 


tions . 


Shell  radius  a 


Distance  between  bulkheads 


Shell  thickness 


5 meters 


16  meters 


0.05  meters 


Circumferential  average  thickness  : 0.07  meters 


Bulkhead  thickness 


Structural  material 


Yield  stress  for  steel 


0.02  meters 


Steel 


700  Megapascals 


Shell  flexural  rigidity  ratio  Dg/D^  : 100 
Shell  twisting  rididity  ratio  Dg^/D^:  0.5 


E.  FAILURE  CRITERIA 

Failure  criteria  used  by  Newton  and  Atchison  are  employed.  The 
following  excerpt  from  reference  [2]  defines  the  basis. 


The  largest  calculated  elastic  stresses  are  found  to  be 
in  the  shell  at  the  bulkhead.  These  stresses  result  from 
superposition  of  longitudinal  bending  amd  axial  membrane 
compression.  Since  yielding  at  this  location  would  not  pre- 
cipitate shell  collapse,  combinations  of  bending  moment  and 
axial  force  corresponding  to  failure  are  based  on  limit 
analysis  (full  yielding  through  thickness  - part  con^ression, 
remainder  tension) . 

A different  situation  exists  at  appreciable  distances 
from  the  bulkhead.  Here  it  is  postulated  that  yielding  will 
induce  shell  collapse.  At  such  locations  em  effective  uni- 
axial stress  a is  calculated  using  the  von  Mises  yield  condi- 
tion. Failure  is  deemed  to  occur  when  the  utilization  ratio 
R =•  a/0  reaches  unity. 


IV.  RESULTS 

A.  INTRODUCTION 

Results  reported  here  concern  two  separate  investigations.  The  first 
of  these  is  a limited  examination  of  the  effectiveness  of  the  radiation 
boundary  condition.  The  second,  and  principal,  study  is  a systematic 
examination  of  the  effects  of  shock  pulse  rise  time  and  duration  on 
structural  failure. 

B.  EFFECTIVENESS  OF  RADIATION  BOUNDARY 

The  effectiveness  of  the  radiation  boundeury  condition  is  examined 
by  repeated  solutions  of  the  same  problem  using  successive  reductions  in 
the  outside  radius  r^  of  the  fluid  region.  Table  I gives  the  correspond- 
ing utilization  ratios  for  both  failure  modes.  The  first  line  of  the 
table  provides  the  comparison  stcuidard, since  the  extreme  values  of  the 
utilization  ratios  occur  at  15.6  milliseconds  cuid  the  time  for  an 
acoustic  wave  to  travel  from  shell  to  the  outer  boundary  and  back  is  18 
milliseconds.  The  second  line  (r^  =>  16.025  meters)  also  represents  a 
ceise  where  the  outer  boundary  is  "far  enough"  from  the  shell  so  that 
the  utilization  ratios  are  unaltered.  Succeeding  radius  reductions 
do  affect  the  utilization  ratios,  but  even  for  the  smallest  region  (r^ 

» 8.15  meters)  the  ratios  chcuige  by  less  than  3 percent.  The  acccanpany- 
ing  reduction  in  the  number  of  fluid  nodes  is  dramatic.  In  reducing 

from  r > 16.025  meters  to  r » 8.15  meters,  the  number  of  fluid  nodes 
o o 

goes  from  120  to  45. 
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TABLE  I.  Effect  of  Radiation  Boundary  Placement  on  Utilization  Ratios 


OUTSIDE  RADIUS 
(METERS) 

TRAVEL  TIME 
(MS) 

RELATIVE  UTILIZATION  RATIOS 

LIMIT  ANALYSIS 

VON  MISES 

17.6 

18.0 

1.000 

1.000 

16.025 

15.8 

1.000 

1.000 

14.45 

13.5 

1.000 

0.998 

12.875 

11.2 

0.997 

0.997 

11.3 

9.0 

0.988 

0.976 

9.725 

6.8 

0.999 

0.974 

8.15 

4.5 

0.995 

0.971 

Notes; 

1.  Shell  radius:  a ^ 5 meters. 

2.  Cubic  ramp-step  wave  with  cT^/a  = 0.97. 

3.  Submergence  pressure  =4.0  Megapascals, 

shock  pressure  = 4.5  Megapascals. 

4.  All  maxima  occur  between  14.3  and  15.7  milliseconds  after  onset. 

5.  Travel  time  is  time  for  acoustic  wave  to  travel  from  shell  to 
outer  radius  and  back  to  shell  (c  = 1.4  meter/millisecond). 

TABLE  II.  Effect  of  Fixed  Outer  Boundary  Placement  on 
Utilization  Ratios 


OUTSIDE  RADIUS 
(METERS) 

TRAVEL  TIME 
(MS) 

RELATIVE  UTILIZATION  RATIOS 

LIMIT  ANALYSIS 

VON  MISES 

17.6 

18.0 

1.000 

1.000 

9.725 

6.8 

0.993 

0.949 

8.15 

4.5 

0.992 

0.947 

Notes;  Notes  for  Table  I apply  here 


T2U3le  II  exhibits  the  effect  of  reducing  the  outer  radius  of  the 
region  without  in^osing  the  radiation  condition  at  that  boundary.  It 
is  appeurent  that  the  resulting  degreuiation  of  accvuracy  is  more  severe 
than  with  the  radiation  condition. 

C.  INFLUENCE  OF  PULSE  RISE  TIME  AND  DURATION  ON  STRUCTURE  FAILURE 

Considered  here  is  the  principal  application  of  the  improved  analy- 
sis techniques  introduced  in  this  thesis.  In  reference  [2]  results 
were  obtained  for  structural  response  to  a cubic  raii5)-step  shock  wave. 

A single  rise  time  cT^/a  ^ 0.97  was  considered  and  the  duration  was 
effectively  infinite.  A large  family  of  cubic  ramp-box  pulses  is 

treated  here.  The  rise  time  T 2md  effective  duration  T » T +T^  + 

r 2 r d 

•j  Tg  eure  systematically  varied.  (Note  that  the  impulse  of  the  shock 
pulse  is  p*T.  See  Fig.  5.) 

Pi 

P* 


Fig.  5.  Effective  pulse  duration  T. 

The  pulse  durations  studied  here  aze  cT/a  « 2,4,  and  8.  For  each 
duration,  four  sepeurate  rise  times:  cT^/a  - 0.1, 0.2, 0.4,  and  0.8  are 
considered.  In  every  instance  the  fluid  region  extends  from  a « 5 
meters  to  r^  ■ 16.025  meters.  For  all  combinations  with  cT^/a  > 0.1  a 
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fluid  mesh  with  90  fluid  elements  and  120  fluid  nodes  is  used  ("coaurse 


mesh).  The  thickness  (radial)  of  each  element  is  0.525  meters.  With 
this  mesh  a time  step  h « 0.3  millisecond  is  used  and  the  subsidence 
time  is  cT^/a  ■ 0.2.  For  the  rise  time  cT ^a  » 0.1  a "fine"  mesh  is 
introduced  with  174  fluid  elements,  225  fluid  nodes  and  a 0.2625  meter 
element  thickness.  With  this  mesh  a time  step  h « 0.1  millisecond  is 
used  and  the  subsidence  time  is  reduced  to  ct  /a  ^ o.l. 

For  both  meshes  there  aure  eight  shell  elements  with  lengths  varying 
from  0.5  meters  near  the  bulkhead  to  2.0  meters  midway  between  bulkheads. 
The  number  of  structural  degrees  of  freedom  is  27  for  each  harmonic. 
Along  the  fluid-structure  interface  there  is  a fluid  node  coinciding 
with  each  structural  node. 

Results  are  presented  graphically  in  the  form  of  failure  boundaries 
on  the  submergence  pressure  versus  shock  pressure  plane  in  Figs.  6 
through  17.  Numerical  results  used  to  define  these  boundaries  are  tabu- 
lated in  Appendix  B. 

Fig.  7 presents  two  sets  of  failure  boundaries  for  cT^/a  « 0.2, 
cT/a  >2.  For  each  boundeury  the  narrow  line  is  found  using  the  fine 
mesh  and  the  wide  line  is  obtained  from  the  coarse  mesh. 

When  the  effective  pulse  duration  is  small  (cT/a  » 2) , changes  in 
pulse  rise  time  affect  considerekbly  the  failure  boundary.  Increase  in 
effective  pulse  duration  (cT/a  ^ 4)  produces  failure  boundaries  that  are 
not  sensitive  to  changes  in  rise  time.  In  fact,  for  cT/a  > 8 the  seuae 
values  are  obtained  for  all  rise  times  greater  than  cT^a  > 0.1.  It  is 
considered  that  the  small  difference  observed  for  cT^/a  - 0.1  is  laurgely 
due  to  the  mesh  change.  For  durations  greater  than  cT/a  > 4 results 
approach  those  of  Fig.  16,  regardless  of  rise  time.  These  results  are 
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based  on  cT^/a  = 0.97  and  cT/a  » <*>,  the  pulse  parameters  used  by  Newton 
and  Atchison  [2] . The  boundaries  found  here  agree  closely  with  those 
of  reference  [2] , but  it  is  noted  that  labels  on  ordinate  £Uid  abscissa 
have  been  mistakenly  interchcuiged  on  Fig.  10  of  [2]. 
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It  is  concluded  that  the  separation  of  the  total  dynamic  pressure 
acting  on  a structure  into  incident  and  scattered  pressures  has  im- 
proved the  analysis  of  the  fluid-structure  interaction  problem  investi- 
gated in  this  thesis  by  avoiding  the  propagation  solution.  The  results 
obtained  also  confirm  that  the  imposition  of  the  radiation  bound^u:y 
condition  at  the  cylindric  outer  boundeury  affords  a further  improvement. 
For  a representative  structure,  the  effects  on  failure  of  pulse  rise 
time  emd  pulse  duration  have  been  established. 
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APPENDIX  A.  FLUID  COEFFICIENT  MATRICES 


1.  Separation  of  fluid  equations 

The  governing  equation  for  the  scattered  pressure  is 

••  • T ** 

2^£+D£+Il£=PL  ~ ^ 


(8) 


where  £ is  the  vector  of  nodal  scattered  pressures.  It  is  assumed 
that  the  (scattered)  pressure  p at  a point  (r,  0,  z)  of  the  fluid  is 
given  by 


P = Nj  Pj  cos  n6 
i in 


(Al) 


where  N^  , a function  of  r and  z , is  the  shape  function  associated 
with  fluid  node  i and  p^^  is  the  nth  Fourier  component  of  pressure  at 
that  node.  Summation  on  indices  i amd  n is  implied. 

This  leads  to  a sepaurate  equation  like  equation  8 for  each  Fourier  com- 
ponent of  vector  £.  The  typical  equation  may  be  written  as 


0£^  + D ^+H  p*pl  6 - f 

-n*Ti  -n^Ti  -n— n “n 


(A2) 


where  the  summation  convention  does  not  apply.  Equation  A2  assumes  that 
Fourier  resolution  has  also  been  applied  to  the  structural  nodal  displace- 
ments and  that  the  normal  displacement  w of  the  shell  is  given  by 


w ■ N ' 6.  cos  n6 
i in 


(A3) 
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where  N^'  , a function  of  z , is  the  shape  function  for  structural 
node  i and  6^^  is  the  nth  Fourier  component  of  displacement  at  this 
node.  As  in  eguaticm  (Al)  the  summation  on  i and  n applies. 


2.  Formulas  for  matrix  elements 

Coefficient  matrices  eure  found  by  calculating  contributions  at  the 
element  level  and  assembling  by  addition.  Presented  in  this  section 
are  the  element- level  contributions  based  on  the  convention  that  a lower 
case  letter  represents  em  element  of  the  matrix  denoted  by  the  corre- 
sponding capital  letter.  Thus,  q . , is  the  element  in  row  i,  column  j, 

ijn 

of  matrix  Q^. 

In  each  of  the  following  formulas  a symbol  e , the  Neumann  factor, 

n 

appears . It  is  defined  to  have  the  value  1 for  n«0  and  the  value  2 for 
n any  positive  integer.  In  the  two  following  formulas  the  integrations 
with  respect  to  r and  z extend  over  the  range  for  the  individual  fluid 
element. 


^ijn 


21T 


e c 
n 


2 //  N^  N^  r dr  dz 


3n.  3N.  3n.  3n.  2 

h . . ■ — fl  (-5 — -5-^  + ^ N N , ) r dr  dz 

ijn  e 3r  3r  3z  3z  2 i j 

n r ■’ 


(A4) 


(AS) 


No  summation  on  n is  intended  in  equation  (AS) . 

The  following  two  formulas  concern  the  boundary  condition  at  the 
structure-fluid  interface  (r«a).  The  only  fluid  nodes  having  non-zero 
contributions  are  those  in  contact  with  the  shell.  Integrations  with 
respect  to  z ta)ce  place  along  a structural  element  and  the  matching 
fluid  element. 
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/ N. ' N . dz 
1 D 


(A6) 


1.  . 

ijn 


b. 

3^ 


n 


N.  dz 


(A?) 


The  component  , when  multiplied  by  as  defined  in  equation  (23) , 

gives  the  jth  component  of  vector  - ^ . 

The  final  two  formulas  represent  contributions  due  to  the  radiation 
boundary  condition  applied  at  the  outer  boundary  (r  * r^) . Integration 
with  respect  to  z is  carried  out  along  that  boundeury.  Only  fluid  nodes 
on  that  boundary  aure  included 


A m 

2iTr 

o 

• / N.  N 

^ j 

t KQ\ 

ijn 

ce 

n 

CLZ 

(Ao; 

Ah.  . * 

ijn 

TT 

e 

n 

/ N,  N. 
i 3 

dz 

(A9) 

The  contributions 

% 

are  assembled  into  H 
n -n 
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APPEtn)IX  B 

NUMERICAL  RESULTS  FOR  FAILURE  BOUNDARIES 


44 


i I 


I ^ 


I i 


j 


FAILURE  BOUNDARY  VALUE 


CTr  - _ 


2 


Number  of  fluid  elements:  174 
LIMIT  ANALYSIS 


shcx:k  pressure 

(MEGAPASCALS) 

SUBMERGENCE  PRES. 
(MEGAPASCALS) 

MAXIMUM  RATIO 

TME  MAXIMUM  RATIO 
(MILLISECONDS) 

10.00 


10.84 


0.999 


1.001 


1.000 


0.999 


1.000 


1.000 


1.000 


SHOCK  PRESSURE 
(MEQAMSCALS) 

SUBMERGENCE  PRES. 
(MEGAPASCALS) 

MAXIMUM  RATIO 

TIME  MAXIMUM  RATIO 
(MILLISECONDS) 

10.20 


8.28 


5.64 


1.84 


1.001 


0.999 


1.001 


1.000 


1.000 
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LIMIT  ANALYSIS 


Number  of  fluid  elements:  i74 

! 

! LIMIT  ANALYSIS 


0.4 


Number  of  fluid  elements:  so 
LIMIT  ANALYSIS 


SHOCK  PRESSURE 
(MEGAPASCALS) 


MAXIMUM  RATIO  TME  MAXIMUM  RATIO 
MAXIMUM  RATIO  (MILLISECONDS) 


0.00 


2.00 


4.00 


6.00 


8.00 


10.00 


10.44 


0.999 


0.999 


1.000 


1.000 


1.000 


1.000 


1.000 


SHOCK  PRESSURE 
(MEGAPASCALS) 


MAXIMUM  RATIO  MAXIMUM  RATIO 

MAXIMUM  RATIO  (MILLISECONDS) 


10.20 

1.001 

8.26 

1.000 

5.63 

1.000 

2.14 

1.000 

0.00 

1.000 
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AILURE  BOUNDARY  VAL 


0.8 


^=:  2 


Number  of  fluid  elements:  9o 
IMIT  ANALYSIS 


SUBMERGENCE  PRES. 
(MEGAPASCALS) 

MAXIMUM  RATIO 

0.00 

9.27 

0.999 

2,00 

7.63 

1.001 

4.00 

5.96 

1.000 

6.00 

4.30 

1.000 

8.00 

2.48 

1.000 

10.00 

0.00 

1.002 

(MILLISECONDS) 


."MOCK  PRESSURE 
\<v£QAPASCALS) 


MAXIMUM  RATIO  TIME  MAXIMUM  RATIO 
MAXIMUM  RATIO  (p*,LLISECONOS) 


0.00 

10.20 

1.001 

• 

O 

o 

8.17 

1.000 

4.00 

5.50 

1.000 

6.00 

1.86 

1.001 

6.76 

0.00 

0.999 
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FAILURE  BOUNDARY  VALUE 


^ = 0.1 


e = 4 


Number  of  fluid  elements: 
LIMIT  ANALYSIS 


FAILURE  BOUNDARY  VALUE 


CTr  _ 


“T  ” °*2 


4L  = 4 


Number  of  fluid  elements:  9o 
LIMIT  ANALYSIS 


SHOCK  PRESSURE 
(MEGAPASCALS) 

SUBMERGENCE  PRES. 
(MEGAPASCALS) 

MAXIMUM  RATIO 

TME  MAXIMUM  RATIO 
(MILLISECONDS) 

0.00 

9.27 

0.999 

2.00 

7.44 

1.000 

4.00 

5.59 

1.000 

6.00 

3.74 

1.000 

8.00 

1.92 

1.003 

8.65 

0.00 

1.001 

SHOCK  PRESSURE 
(MEGAPASCALS) 

SUBMERGENCE  PRES. 
(MEQAMSCALS) 

MAXIMUM  RATIO 

0.00 

10.20 

1.001 

2.00 

7.90 

1.001 

4.00 

4.90 

1.000 

6.00 

0.90 

1.000 

6.35 

0.00 

1.000 

(MILLISECONDS) 


FAILURE  BOUNDARY  VALUE 


0.4 

CT  « 4 

a ” 


Number  of  fluid  elements:  9o 
LIMIT  ANALYSIS 


SHOCK  PRESSURE 
(MEGAPASCALS) 

SUBMERGENCE  PRES. 
(ME(3APASCALS) 

MAXIMUM  RATIO 

TIME  MAXIMUM  RATIO 
(MILLISECONDS) 

0.00 

9.27 

2.00 

7.44 

4.00 

5.60 

6.00 

3.76 

8.00 

1.92 

8.69 

0.00 

0.999 


1.000 


1.000 


1.000 


1.001 


13.2, 

14.4 

13.2, 

14.4 

13.2, 

14.4 

13.2, 

14.4 

13.2, 

14.4 

SUBMERGENCE  PRES. 
(MEGAPASCALS) 

MAXIMUM  RATIO 

10.20 

1.000 

7.90 

1.000 

4.94 

1.001 

0.90 

1.001 

0.00 

1.000 

(MILLISECONDS) 


13.2,  14.4 


14.4 


14.4 


14.4 


Number  of  fluid  elements: 


von-MIses 


Number  of  fluid  elements:  i74 


LIMIT  ANALYSIS 


SHOCK  PRESSURE 
(MEGAPASCALS) 

SUBMERGENCE  PRES. 
(MEGAPASCALS) 

MAXIMUM  RATIO 

0.00 

9.27 

0.999 

2.00 

7.37 

1.000 

4.00 

5.46 

1.000 

6.00 

3.55 

1.000 

8.00 

1.26 

1.000 

8.38 

0.00 

1.000 

(MILLISECONDS) 


SHOCK  PRESSURE 
(MEGAPASCALS) 

SUBMERGENCE  PRES. 
(MEGAMSCALS) 

MAXIMUM  RATIO 

0.00 

10.20 

1.001 

2.00 

7.81 

1.001 

4.00 

4.78 

1.000 

6.00 

0.78 

1.003 

6.08 

0.00 

1.000 

(MILLISECONDS) 


12.6 


12.6,  13.2 


13.2 


13.2 
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Number  of  fluid  elements: 


FAILURE  BOUNDARY  VALUES 

^ = 0.97 

^ = ~ 

Number  of  fluid  elements.  ^02 


LIMIT  ANALYSIS 


SHOCK  PRESSURE 

SUBMERGENCE  PRES. 

MAXIMUM  RATIO 

TME  MAXIMUM  RATIO 

(MEGAPASCALS) 

(ME(3APASCALS) 

(MILLISECONDS) 

0.00 

9.27 

0.999 

2.00 

7.45 

1.000 

14.4,  15.6 

4.00 

5.62 

1.000 

14.4,  15.6 

6.00 

3.75 

1.000 

14.4 

8.00 

1.95 

1.000 

14.4 

8.70 

0.00 

0.998 

14.4 

von-MIses 


SHOCK  PRESSURE 

MAXIMUM  RATIO 

TIME  MAXIMUM  RATIO 

(MEGAMSCALS) 

■rv  1 iVi' 

(MILLISECONDS) 

0.00 

10.20 

1.001 

2.00 

7.91 

1.000 

14.4 

4.00 

4.95 

1.000 

14.4 

6.00 

0.94 

1.000 

15.6 

S6 
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